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1. Introduction 

In recent years, precise measurements of tlie Cosmic Microwave Background (CMB), 
large scale structure (LSS), and distant supernovae (SNe) have increased our knowledge 
of the large scale properties of the universe tremendously. The standard ACDM model of 
cosmology has turned out to be able to fit all available data with only a few free parameters. 

However, there are still many parameters related to cosmology which are not well de- 
termined by the present large scale observations. For instance, in order for the CMB data 
from the WMAP experiment (Spergel et al. 2003; Bennett et al. 2003; Kogut et al. 2003; 
Hinshaw et al. 2003; Verde et al. 2003; Peiris et al. 2003) to provide a stringent limit on 
the curvature of the universe, a prior on the Hubble parameter is necessary. This degeneracy 
can be broken by adding LSS data from surveys such as the Sloan Digital Sky Survey (SDSS) 
(Tegmark et al. 2004a; Tegmark et al. 2004b) or the 2dF survey (Peacock et al. 2001), but 
other independent methods for measuring the Hubble parameter precisely are of great inter- 
est. 

Another problem is that the nature of the dark matter is at present unknown. The 
properties of dark matter on large scales are compatible with a heavy, coUisionless species 
which clusters gravitationally, cold dark matter (CDM). On the other hand there are signif- 
icant problems with understanding some properties of galaxies in CDM models. Numerical 
simulations of galaxy formation in CDM models show that dark matter halos exhibit a uni- 
versal density profile with a central slope p oc r"'', with r] ~ 1 — 1.5 (Navarro et al. 1996; 
Ghigna et al. 2000; Power et al. 2003; Fukushige et al. 2004; Hayashi et al. 2003; Kazantzidis et al. 2004; 
Navarro et al. 2004). At large radii the profile steepens to p oc r~^. While the outer slope is 
consistent with current observations, there are problems with fitting the inner slope of halos. 

Observations of rotation curves of many dwarf galaxies, indicating almost constant den- 
sity cores, suggest that the inner density profiles in these systems are much shallower than 
found in simulations (de Blok et al. 2001; van den Bosch et al. 2000; Dutton et al. 2003; 
Simon et al. 2003a; Simon et al. 2003b). Another problem is that the number of subha- 
los within a Milky Way-sized halo predicted by the simulations exceed the observed number 
of satellite galaxies by at least an order of magnitude. 

Many explanations have been offered for this discrepancy. One possibility is that tidal 
interactions between baryons and dark matter erase the central cusps, and that early star 
formation expels baryons from small subhalos, rendering them invisible. Whether this purely 
astrophysical explanation can work is not clear at present. Another possibility is that the 
properties of dark matter are slightly different from pure CDM. If for instance dark matter 
has relatively strong self-interactions or is warm instead of cold, shallow density profiles and 



-3- 



lack of subhalos might be explained. However, it should be noted that the simplest models 
of self-interacting dark matter predict large constant-density cores in clusters of galaxies 
which are ruled out by X-ray and gravitational lensing observations (Meneghetti et al. 2001; 
Arabadjis, Bautz, & Garmire 2002; Dahle, Hannestad, & Sommer-Larsen 2003). 

At present the ultra-high resolution simulations of galaxy formation do not include 
baryons. This in turn means that for regions of galaxies which are baryon dominated the 
discrepancy might be resolved by a more complete understanding of feed-back from star 
formation. In dwarf galaxies this effect is unlikely to be dominant, but for the high-mass 
systems typically investigated with the method proposed here, the regions probed will be 
dominated by the baryonic component. Therefore the density profile which we derive for 
the inner parts of galaxies is likely to be dominated by baryonic physics and can therefore 
not be directly related to either the pure N-body simulations or to the observations of dark 
matter dominated dwarf galaxies. 

Nevertheless it is of great importance to get a precise observational determination of 
the halo density profiles in order to understand whether the difference between observations 
and simulations is generic, or specific to some systems. 

Gravitational lensing has the potential to give information on both the matter distribu- 
tion in galaxies as well as the Hubble parameter. More than 40 years ago, Refsdal outlined 
how it should be possible to measure the Hubble parameter and the mass of a galaxy by 
measuring the time delay between multiple images of a supernova (SN) lying far behind and 
close to the line of sight through the lensing galaxy (Refsdal 1964). 

The first observation of gravitational lensing in a cosmological context was made in 1979, 
when two images of the quasar QSO 0957+651 at 2; ~ 1.4 was observed. For several reasons, 
it has been difficult to implement Refsdal's method on this system. Until 1997 it was not 
possible to precisely determine a time delay for the system, and even then the complexity 
of the lens system (a brightest cluster galaxy located close to the center of a galaxy cluster) 
has made estimates of the Hubble parameter highly uncertain (Keeton et al. 2000). Studies 
of this system have also showed the importance of including effects from substructure in the 
lens. 

To date, approximately 70 multiply imaged sources have been observed out of which 
there are eleven firm time delay measurements (Davis et al. 2003; Kochanek & Schechter 2003; 
York et al. 2004). There is not yet any agreement on the correct value of the Hubble pa- 
rameter inferred from these measurements due to differences in the modeling of the lensing 
galaxies. In general, modeling the lenses as more concentrated gives larger values of h and 
vice versa, e.g., Kochanek (2002) obtains h ~ 0.48loo4 if the lenses have isothermal mass 
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distributions and h ~ 0.71 ± 0.06 if they have constant mass-to-hght ratio. 

There are basically two different routes in trying to improve on current results. One is to 
try to do extremely careful modeling of single simple lens systems. Following this approach, 
Wucknitz, Biggs, & Browne (2004) find h = 0.78 ± 0.06 (2 a) based on the lens system 
B0218+357, consistent with the local estimate h = 0.71 ± 0.06 (1 a) from the Hubble Space 
Telescope (HST) key project (Mould et al. 2000; Freedman et al. 2001). Another possibility 
is to combine results from a larger number of, perhaps less well constrained, systems. In this 
paper, we will investigate the latter method. 

In order to avoid any systematic bias due to selection effects, it is important to have 
a well-defined statistical sample. The currently largest survey of strong lensing events is 
the recently completed Cosmic Lens All-Sky Survey. A subsample of 8958 fiat-spectrum 
radio sources of which 13 have multiple images constitutes a well-defined statistical sam- 
ple and has been used to constrain cosmological parameters as well as galaxy properties 
(Davis et al. 2003; Chae et al. 2004). 

Future surveys that will improve statistics and also include multiply imaged SNe for 
which time delay measurements can be obtained at high accuracy due to the characteristic 
light curves include Pan-STARRS ^ and the Supernova Acceleration Probe (SNAP) ^. 

In this paper we investigate how future SNAP data of multiply imaged core-collapse 
SNe can be used to measure the slope of galactic halos and the Hubble parameter (see also 
Holz 2001; Oguri et al. 2003). This paper extends and generalizes earlier work 
(Goobar et al. 2002a) where similar data was used to constrain h and Qu using simpler 
lens models. A recent work by Oguri & Kawano (2003) have considered a similar use of 
quadruply lensed Type la SNe events. The possiblity of using SNe that are multiply lensed 
by rich galaxy clusters to constrain h has been considered by Bolton & Buries (2003). The 
complementarity of strong lensing to other probes of the expansion rate is discussed by 
hinder (2004). 

Finally, it should be noted that while the method proposed here is very sensitive to 
the inner density profile of galaxies, this density profile does not necessarily refiect the 
nature of pure dark matter halos. Lensing mainly occurs in relatively massive galaxies in 
which baryonic material is a dominant component at small radii. This in turn means that 
there can be a substantial modification of the inner density profile due to cooling and infall 
(Kochanek & White 2001). 



^http : //www . pan-Starrs . org 

'"'SNAP Science Proposal, available at http://snap.lbl.gov 
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In Sec. 2, we describe the lens model used in the simulations and analysis and in Sec. 3 
we describe our numerical simulations. Errors are discussed in Sec. 4 and our results are 
presented in Sec. 5. Finally, Sec. 6 contains a summary of our results. 



2. Method 

2.1. Lens model 

For most lensing purposes, galaxy matter distributions can be described by the projected 
Newtonian gravitational potential 

We assume that lens systems are characterized by a simple power-law density profile 

p{r) oc r-". (2) 

This is not a good approximation at all radii since the slope is expected to change with 
radius. In practice however, almost all images of multiply lensed SNe are at a limited range 
of small r (see Fig. 1) and thus only the inner slope is probed and the single power-law 
density profile is an excellent approximation. 

For the power-law density profile, the lenses can be described by the projected potential 

^ = l^x'-^ (3) 

3 — 1] 

where x is the impact parameter (in arbitrary units, .^o) and xe is the Einstein radius. 



2.2. Time delay 

The time delay for a gravitationally lensed image as compared to an undefiected image 
is in the general case given by (Schneider et al. 1992) 






y)' 



(4) 



where ^o is the (arbitrary) scale length, the D's are angular diameter distances and y is the 
source position. The time delay between two lensed images is given by At = T2 — Ti. For 
an isothermal lens, we have rj = 2 and 

Ats,s = §^a + zi){xl-xl). (5) 
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Putting ^0 = Di, i.e., denoting positions in terms of angles we get 

IDiDs 



At 



SIS 



2 D 
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'■-{i + z,){ei-e',). 



(6) 



In the general case of rj ^ 2, we can rewrite Eq. (4) in terms of g = A6/ <6> where A^ = 
9i — 02 and <9>= {9i + ^2)/2 and Taylor expand to get [after some tedious calculations; see 
also Kochanek & Schechter (2003) and Chang & Refsdal (1977) for an earlier approximation] 
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(7) 



We see that we have an almost perfect degeneracy between rj — 1 and h (which comes in 
through the distances). This degeneracy can be broken by including data from the observed 
flux ratio. 



2.3. Flux ratio 

We denote the flux ratio r = /ii/|/i2|. The magnification fi for a source at position y 
observed at position x is given by 



l/^l 



ydy 



Computing the lensing angles through a = V^ and using the lens equation 



we can express the flux ratio as 



X2 

^sis 



y = xi - ai = X2- a2, 



l-{2-r])x2 '^{Xi + X2)/{X;^ "^ + X2 '') 



2-r] 
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l-{2-r])xi '^{xi+ X2)/{x^ "^ + X2 '^) 
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l-{2-r])x^ '^{xi + X2)/{x^ '^ + X2 '') 

Expanding this to second order in q, we obtain 

1 
rc^l + {ri-l)q+-{ri- 1) V + 0{q^). 



(9) 



(10) 



;ii) 



Eqns. (7) and (11) are excellent approximations for small q and/or r^ ~ 2. In order 
to show the general parameter dependencies, we generalize in the following the expanded 



-7- 



expressions to include also the effects from ellipticities and external shear. However, in the 
subsequent x^-analysis, we use the full expressions for the time delay and the flux ratio to 
avoid any bias due to the Taylor expansion^. 



2.4. Ellipticity 

In order to treat possible non-sphericity of the lens systems we add a quadrupole term 
to the projected potential, \E', 

^ = ^0(1 + a cos 20), (12) 

where \E'o is the spherical potential given in Eq. (3) and is the angle of the image relative 
to the quadrupole axis. Using this potential it is possible to derive an analytic expression 
for the time delay, which is, however, quite complicated. If the system is assumed to have 
relatively small ellipticity, an assumption already implicit in the fact that we represent the 
asphericity of the lens with a quadruple moment, then the angle of the second image can be 
written as a function of the first 

02 = 0i + vr-5, (13) 

where (5 is a small parameter (and explicitly for spherical systems). We expand the time 
delay in the parameters q and 5 to find 



Atc^{ri- l)At; 



SIS 



12 ^ 4 



(14) 



an expression which is valid to second order in both q and 5. In the same way it is possible 
to derive an expression for the flux ratio of a system with a quadrupole moment. To second 
order in q and 5 we find that 

3 + cos(40i) ^r, 
r = l + 2^cot(20O+ . 2,0.2? ^ 

sm (20f) 

HV - 1) [1 + 5cot(20i)] q+\{v- 1)'?' + 0{q\ 5\ q\ q5'). (15) 

The above expression appears divergent for 0i -^ because of the l/sin(20i) terms. How- 
ever, since 5^0 when 0i -^ 0, the quantity (5/sin(20i) remains finite, and there are no 
actual divergences. 



^Since 77 ^ 2 in our simulations, the use of the Taylor expanded expressions gives identical results. 



2.5. External shear 



Most lens systems are embedded within an external potential which gives rise to an ad- 
ditional shear contribution. Following Kochanek & Schechter (2003) we write the combined 
potential as 

^ = ^0 (1 + a cos 20) + 7extr'cos[2(0 - 0o)], (16) 

where — 0o is the angle between the image position and the quadrupole axis of the external 
potential. If the external shear is aligned with the lens system itself then 0o = 0, but in the 
general case 0o 7^ 0. 

Since 7 should be small in order to justify treating the external potential as a quadrupole 

(an assumption supported in, e.g., Dalai & Watson 2004), we derive expressions for At and 

r which are valid to second order in g, 6 and 7 for arbitrary values of 0o- The expressions 

we find are 

{v-2? . . 1 



and 
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sin(0i)2 cos(0i) 



^^, ^, 16(5(r7-l)^cos(20o)sin(20i 
lb{ri — 1 

2 ^ -6[81 + {r]- 98)r]] cos(20i) + (?] - 1) [6{r] - 1) cos(60i) + 8qr] sin(20i)] 



3. Simulated data 

With the proposed SNAP satellite a large number of core collapse supernovae (CC SNe) 
will be discovered, of which a small fraction will be multiply imaged. We use the SNOC 
package (Goobar et al. 2002c) to simulate a total of 1.1 • 10^ CC SNe, the predicted number 
following the prescriptions in Dahlen & Fransson (1999) for a 20 square degree field during 
three years for z < 5. Lensing effects are investigated by ray-tracing using SIS halo profiles 
since lensing statistics are insensitive to ellipticities and shear (Huterer 2004). We derive 
a galaxy mass function by combining the Schechter luminosity function with a mass-to- 
luminosity relation for fundamental plane ellipticals; see Bergstrom et al. (2000) for details. 
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The cosmology used is h = 0.65, ^m = 0.3 and ^a = 0.7. Out of the 1.1 ■ 10^ CC 
SNe, 2613 have multiple images [this is the sample used in Goobar et al. (2002a)]. Out 
of these, 857 have either an I-band or J-band peak brightness < 28.5 mag for the dimmest 
image. In order to avoid contamination from the lens galaxy in the SN detection, we impose a 
constraint that the surface brightness of the lens galaxy has to be fainter than 24 I magnitudes 
per square arcsecond. We assume that the lens galaxies follow a de Vaucouleurs profile 
(de Vaucouleurs 1948) 

lir) = h exp {-7.67 [{r/r^f^ - l] } (19) 

where I{r) is the surface brightness at radius r, r^ is the half light radius and Jg is the surface 
brightness at r,,. Observationally, the surface brightness at r,, (in B magnitudes per square 
arcsecond), /xse, can be related to the half light radius according to (Kormendy 1977) 

/iBe = 3.02 log(re/l kpc) + 19.74. (20) 

The absolute B-band magnitude, Mb, can be written 

Mb = /iBe - 5 log(re/l kpc) - 39.28. (21) 

We use Eqns. (19)-(21) together with the Faber- Jackson relation L oc f ^ and a K-correction 
appropriate for early-type galaxies to calculate the I-band surface brightness, /ij at the posi- 
tion of the faintest image. Here, we have neglected the modest surface brightness evolution 
of elliptical galaxies to 2; ~ 1. All but 4 of the 857 SNe satisfy /ii > 24 at this position and 
we conclude that the lens galaxy light does not significantly affect the SN observations. 

Since we assume fairly clean and simple lens systems in our analysis, we incorporate a 
quality factor, /, giving the fraction of the lens systems that fulfill our requirements of a 
single dominant deflector with not too large ellipticity or external shear. In accordance with 
currently observed systems (see, e.g., Kochanek & Schechter 2003), we set / = 0.5 in our 
analysis. We also perform an extremely conservative analysis setting / = 0.05. 

Even though we can derive fair lensing statistics using a simple SIS galaxy profile, it is 
important to include lens ellipticities and galaxy environments in the analysis 
(Keeton & Zabludoff 2004). We modify our simple simulated lensing systems by adding the 
effects of ellipticities and shear according to the distributions described in Sec. 4. Since we 
expect a scatter in the properties of individual halos, we have also added a 1 sigma dispersion 
in the value of 77 for individual halos of Ar^ = 0.2. However, in order to constrain the mean 
slope of galactic halos, we are still fitting r] as a global parameter. 
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Assumed errors 



ar'. Errors in the flux ratio of tlie two SN images will be dominated by the effects of 
microlensing by compact objects (Schneider & Wagoner 1987) and millilensing by cold dark 
matter (CDM) sub-halos (Dalai & Kochanek 2002; Keeton 2003; Mao et al. 2004), rather 
than flux measurement uncertainties. 

In QSO microlensing, a source of essentially constant size is observed moving across a 
microlens ampliflcation pattern at a typical transverse velocity of ~ 500 km s~ . In the case 
of microlensing of a SN, this transverse velocity is dwarfed by the ~ lO^kms"^ expansion 
velocity of the SN photosphere. If the center of the SN is close to a critical line, an amplifl- 
cation of several magnitudes may result from microlensing (Schneider & Wagoner 1987). In 
this case, the shape of the light curve of the affected SN image is usually strongly modifled, 
and such cases should be readily identiflable from a comparison of the light curves of the 
two SN images. 

The light curves of QSO 2237+0305 (Huchra's lens) presented by the OGLE team 
(Wozniak et al. 2000) indicate a scatter 0"^ = 0.4 caused by stellar microlensing. When also 
taking into account a ~ 20% uncertainty in the image fluxes due to CDM substructure, we 
arrive at a conservative estimate of (T^/r = 0.5, which was used for our simulations. 

a^: Following Goobar et al. (2002a), we conservatively estimate a positional uncertainty 
of (jQ = 0.01" for SNAP data. While millilensing by CDM sub-halos may affect the flux 
ratio, its effect on image positions would likely be negligible, compared to the measurement 
uncertainty. 

At: The estimated measurement uncertainty for the time delay is 0.05 days (Goobar et 
al. 2002a). However, in the presence of microlensing, the measured time delay may depart 
signiflcantly from the time delay predicted by our smooth macro-lens model. In the most 
extreme cases, a high-ampliflcation microlensing event may shift the observed maximum by 
several days (Schneider & Wagoner 1987), but as noted above, such rare events should be 
identiflable and can thus be removed from the sample of events. We model the observed 
SN light curve as the intrinsic curve superposed on a smooth gradient, with a typical value 
similar to the microlens light curves of Wozniak et al. (2000). From this, we estimate an 
uncertainty a At of 0.15 days. 

7cxt: Using N-body simulations and semi-analytic models of galaxy formation. Holder & 
Schechter (2003) estimate a mean external shear at the position of lens galaxies of 7ext = 
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0.058, with an rms dispersion of 0.071. They also find a significant tail towards high values 
of 7cxt- We note that the photometric redshift estimates obtained with SNAP should enable 
the identification - and removal - of a small subset of lens systems where nearby galaxy 
groups or clusters produce large values for the external shear. 

Quadrupole: Detailed models of known gravitational lens systems indicate that the mass 
distribution is aligned with the optical light distribution, with an upper limit (A0^)^/^ < 10° 
on the dispersion in the angle between the major axes (Kochanek 2002). The measurement 
uncertainties in determining the positions of the lensed images with respect to the lens galaxy 
will be negligible by comparison. In lens galaxies undergoing major mergers, a temporary 
misalignment of the optical and dark matter may occur (Quadri, Moller & Natarajan 2003). 
Here, we neglect such effects on the assumption that such merging systems will have mor- 
phological and photometric signatures that enable their identification and removal from our 
sample. 

azi'. We adopt o-^^ = 0.001 from photometric redshift measurements based on SNAP multi- 
band photometry (Goobar et al. 2002a). 

fi^: We use fi„ = 0.30±0.04 from the Tegmark et al. (2003b) analysis of SDSS and WMAP 
data. It should be noted that this is a very conservative estimate of the error. Future CMB 
measurements from experiments like Planck, as well as large scale weak lensing surveys are 
likely to increase the precision with which Vim can be measured. 



5. Results 

In order to investigate constraints from the time delay measurements, we perform a 
X^-analysis over a [h, ?7]-grid 

X\K V) = [Atr - Atf (/., v)]V,-'[Atf^ - Atfih, v)], (22) 

where Atf^^ is the experimental time delays (simulated using the SIS profile and modified to 
include the effects from ellipticity, external shear and a scatter in the slope, see Sec. 3) and 
Atf^{h,ri) is the theoretical time delays. Vij is the error (covariance) matrix. We expect a 
close to perfect degeneracy between h and l — rj that can be broken by including information 
from the fiux ratio that is sensitive to the slope t] only. For constraints using the fiux ratio, 
we use 

X'ih) = [rr - rnv)n'[rr " ^?(v)], (23) 
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where r^^^ is the experimental (simulated) flux ratio and rf^{i]) is the theoretical flux ratio. 

In order to avoid any bias due to the second order expansion in q in the time delay 
and flux ratio, we use the full expression for these quantities. However, since r^ ~ 2 in our 
simulations, the use of Eqns. (17) and (18) for the theoretical time delay and flux ratio, 
respectively, gives close to identical results. 

In Fig. 4, results from the time delay (left panel), flux ratio (middle panel) and the 
combined results (right panel) are shown. Using a quality factor / = 0.5, this gives a total of 
~ 400 SNe for the magnitude cuts described in Sec. 3. Contours correspond to 68.3 %, 90 %, 
95 % and 99 % confidence levels. It is clear that when combining the results from the time 
delay and flux ratio measurements, we are able to make a determination of h within 10 % 
and, perhaps more interestingly, to determine rj at the per cent level at 95 % confidence. 

In Fig. 5, results using a quality factor / = 0.05, giving a total of ~ 40 SNe are shown. 
Even with this drastic decrease in statistics, we are still able to determine the slope rj to an 
impressive accuracy. 

Constraints from the time delay are almost fully determined by the quality of the time 
delay measurements whereas the fiux ratio constraints are quite robust to all observational 
errors. Even when increasing the error in the fiux ratio observation with a factor of 3 (to 
150%), we are still able obtain 77 = 2 ± 0.06 at 90% confidence. Increasing the size of 
the external shear by a factor of ten causes a systematic bias in the determination of r] of 
?7 = 2 — > 1.65. Increasing the value of the dispersion in the value of 77 for individual halos to 
At] = 0.5 will cause a bias in h oi h = 0.65 — > 0.55 but the mean slope r] will still be well- 
determined. Thus, we conclude that our results are fairly robust even to quite significant 
changes in the quality and quantity of the data used in the analysis. 

As is evident from Eq. (18), the effect on the fiux ratio from varying rj is strongest for 
high q, i.e, large fiux ratios. For large q, we also have large time delays. Since these systems 
have smaller fractional time delay errors, we obtain stronger constraints also from the time 
delay analysis for large q. Thus, we conclude that systems with large fiux ratios and time 
delays (i.e., large q) are very important when constraining galaxy density profiles and the 
Hubble parameter using strong gravitational lensing. This is confirmed by dividing the lens 
system into two classes with < g < 1 and 1 < g < 2 where (for an equal number of 
lens systems) the class with high q gives confidence contours a factor of ~ 2 smaller in the 
/i-direction and a factor of > 5 smaller in the ?7-direction compared with the low-q systems. 

In practice, it should be possible to obtain additional constraints for a significant fraction 
of the observed lens systems, since some of the SN host galaxies will be lensed into multiple, 
resolved images that may be detectable in deep combined frames from the survey. We 
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make no attempt here to model the possible improvements in the parameter determination 
resulting from such additional constraints, but note that our results here correspond to a 
pessimistic scenario where no multiply lensed SN host galaxies are detectable. 



6. Summary 

We have analyzed how the measurement of strongly lensed SNe by future SN surveys 
will allow for a very precise determination of statistical properties of the matter distribution 
in galaxies, as well as the Hubble parameter. 

Understanding the properties of dark matter in details is one of the most important 
problems in modern cosmology, and one of the best laboratories for this is galactic halos. 
While CMB and LSS data probe properties of dark matter on very large scales, measure- 
ments of galactic halos probe small scale properties of dark matter, such as free-streaming, 
dark matter self-interactions (Spergel & Steinhardt 2000; Dave et al. 2001; Hannestad 1999; 
Wandelt et al. 2001; Colin et al. 2002; Burkert 2000; Firmani et al. 2000; Yoshida et al. 2000), 
dark matter - baryon interactions (Chen et al. 2002; Boehm et al. 2002), or modified pri- 
mordial power spectra (Kamionkowski & Liddle 2000; Sigurdsson & Kamionkowski 2004). 

While the method proposed here is not directly applicable to dwarf galaxies where dark 
matter properties can be probed directly it allows for a very precise determination of the 
density profiles of massive galaxies. This in turn can shed light on the feedback between 
baryons and dark matter during the epoch of non-linear structure formation. 

We estimate that with future data it will be possible to constrain the inner density 
profiles of galactic halos to an accuracy of 1-2%, much better than any present measurement. 
Furthermore our proposed method mainly probes lens systems at relatively high redshift, 
z > 1, where no other reliable method for measuring galaxy density profiles is available. 

The estimated precision is based on a sample of roughly 850 multiply imaged SNe. Even 
with a significantly smaller set of SNe, tight constraints could be obtained. This means that 
properties like redshift evolution of density profiles, variations in density profiles with galaxy 
luminosity, color and many other important parameters can be measured by dividing the 
lens systems into different categories. This in turn will provide valuable information on the 
physics of galaxy formation. 
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Fig. 1. — Distribution of impact parameters for primary images for the lens systems used in 
the analysis in this paper. 
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Fig. 2. — Le/t panel: The number of observed multiply imaged SNe as a function of magnitude 
threshold for the dimmest image in the I-band or J-band. Right panel: The solid line shows 
the distribution of source redshifts, the dotted line the lens redshifts for the sample of 853 
SNe that have either an I-band or J-band peak brightness < 28.5 mag for the dimmest image 
and host galaxy surface brightness fii > 24 mag at the position of this image. 
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Fig. 3. — Upper left: The time delay distribution for 853 SNe that have either an I-band or 
J-band peak brightness < 28.5 mag for the dimmest image and host galaxy surface brightness 
fii > 24 mag at this position. Upper right: The distribution of flux ratios for the same sample. 
Lower left: The distribution of g-values. Lower right: The image separation distribution. 
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Fig. 4. — Results from the time delay (left panel), flux ratio (middle panel) and the combined 
results (right panel) using a quality factor / = 0.5 or ~ 400 SNe. Contours correspond to 
68.3%, 90%, 95% and 99% confidence levels. 
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Fig. 5. — Results from the time delay (left panel), flux ratio (middle panel) and the combined 
results (right panel) using a quality factor / = 0.05 or ~ 40 SNe. Contours correspond to 
68.3%, 90%, 95% and 99% confidence levels. Note different axis range from Fig. 4. 



